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Introduction 


The main peptidase (MP"°’) from the coronavirus (CoV) causing severe 
acute respiratory syndrome (SARS) is one of the most attractive molecular 
targets for the development of anti-SARS agents. We report the irreversible 
inhibition of SARS-CoV MP"® by an aza-peptide epoxide (APE; Kinact/ Ki= 
1900(+400) M~'s~'). The crystal structures of the M?’°:APE complex in 
the space groups C2 and P22;2, revealed the formation of a covalent bond 
between the catalytic Cys145 S” atom of the peptidase and the epoxide C3 
atom of the inhibitor, substantiating the mode of action of this class of 
cysteine-peptidase inhibitors. The aza-peptide component of APE binds in 
the substrate-binding regions of M?"® in a substrate-like manner, with 
excellent structural and chemical complementarity. In addition, the 
crystal structure of unbound MP"® in the space group C2 revealed that 
the “N-fingers” (N-terminal residues 1 to 7) of both protomers of MP" are 
well defined and the substrate-binding regions of both protomers are in the 
catalytically competent conformation at the crystallization pH of 6.5, 
contrary to the previously determined crystal structures of unbound MP" 
in the space group P2). 

© 2005 Elsevier Ltd. All rights reserved. 
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and drugs are in progress, these agents are still far 
from clinical use.. Additional efforts in these areas 


Severe acute respiratory syndrome (SARS) first 
emerged in China in November 2002. This highly 
transmissible, infectious and often fatal disease 
spread to 32 countries across five continents, 
causing close to 8500 infections and over 900 deaths, 
until being contained by the summer of 2003. 
Several infections in Asia were reported sub- 
sequently, alerting the world that it remains at risk 
of another outbreak of SARS (World Health 
Organization: Severe acute respiratory syndrome’). 
Although the development of anti-SARS vaccines 


Abbreviations used: SARS, severe acute respiratory 
syndrome; CoV, coronavirus; APE, aza-peptide epoxide; 
MP", the main peptidase from CoV; CMK, chloromethyl 
ketone. 
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of study therefore remain paramount. 

SARS is caused by a novel coronavirus (CoV);** 
it is an enveloped positive-sense single-stranded 
RNA virus infecting respiratory and gastrointesti- 
nal epithelial cells, macrophages, and other cell 
types, thereby causing systemic changes and 
damaging many vital organs such as lung, heart, 
liver, kidney and adrenal gland.*” Anti-SARS 
therapeutics could target several major steps in 
the viral life-cycle, such as virus—cell interactions, 
virus entry, intracellular viral replication, and virus 
assembly and exit.’ The intracellular replication of 
SARS-CoV is mediated by a “replicase” complex 
derived from two virally coded polyprotein_pre- 
cursors, ppla (486 kDa) and pp1lab (790 kDa).°” The 
formation of this replicase complex requires the 
extensive processing of the two polyproteins by two 
cysteine peptidases within them, namely the main 
peptidase (MP), also known as the 3C-like 
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Figure 1. Inhibition of SARS-CoV MP"’ by aza-peptide 
epoxides (APEs). (a) APEs synthesized for our study, 
Cbz-Leu-Phe-AGIn-EP-COOEt. The epoxide carbon 
atoms are numbered and their stereochemistries are 
omitted for simplicity. The proposed mechanism for the 
irreversible inhibition of clan CD cysteine peptidases by 
APEs is indicated by arrows. Cbz, the benzyloxycarbonyl 
eroup; AGIn, aza-glutamine; EP, epoxide; COOEt, ethyl 
ester. (b) Progress curves for the steady-state cleavage of a 
fluorogenic peptide substrate observed using 64 uM 
peptide and either no APE (dotted line) or 5 uM Cbz- 
Leu-Phe-AGIn-(S,S)EP-COOEt (continuous line). The 
appearance of product was followed using excitation 
and emission wavelengths of 320 and 420 nm, respect- 
ively. Analysis of these data using equation (1), in which 
P; and P. represent the initial and final product 
concentrations, _ respectively, yielded inactivation rates, 
js, of 0.07 min~' and 0.29 min ", tespectively. (c) The rate 
of MP*° inhibition was determined using 16 uM (LJ), 
32 uM (Mi), 64 uM (A) and 100 uM (CO) peptide. Equation 
(2) was fit to the data using the least-squares, dynamic 
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peptidase (3CL?"") because of its similarity to the 3C 
peptidases of Picornaviridae,? and the accessory 
papain-like peptidase 2 (PL2P"°). PL2P"° cleaves at 
three sites in the N-proximal regions of the two 
polyproteins, whereas MP™ cleaves at 11 sites in the 
central and C-proximal regions of the two poly- 
proteins. MP" releases the key proteins in viral 
replication, such as the RNA- -dependent RNA 
polymerase and the helicase.’ Playing such an 
essential role, SARS-CoV MP"® is an attractive 
molecular target for the development of anti-SARS 
drugs acting as the inhibitors of the peptidase. 

SARS-CoV MP" has a molecular mass of 33.8 kDa 
per protomer; it exists as a homodimer over a wide 
range of concentrations in solution.”"'* The crystal 
structures of MP" in the space group P2, showed 
that the two protomers of the dimeric peptidase are 
oriented almost perpendicular to each other and 
that each protomer consists of three domains. 
Domain I (residues 8 to 101) and domain II 
(residues 102 to 184) comprise a two-f-barrel fold 
similar to that of the chymotrypsin-type serine 
peptidases. Domain III (residues 201 to 300) has five 
a-helices and is connected to domain II by a long 
loop (residues 185 to 200). Each protomer has its 
own. substrate-binding region situated in the cleft 
between domains I and II.'° A recent mutagenesis 
study has confirmed that MP" is a cysteine 
peptidase with a Cys-His catalytic dyad at the 
active site.'* As suggested by the structure-based 
sequence alignment of the main peptidases 
(including their flanking residues in the polypro- 
teins) from SARS-CoV and other coronaviruses, 
and confirmed by in vitro studies,” these pepti- 
dases preferentially cleave at a consensus sequence 
for the P4 to P1’ residues of substrates (nomen- 
clature based on that of Schechter and Berger’’): 
(amino acid with a small side-chain)-(any amino 
acid)-Leu-Gln | (Ala, Ser, Gly). 

A number of small-molecule inhibitors of SARS- 
CoV MP" have been proposed using various 
methodologies, such as knowledge-based dis- 
covery” ”® and high-throughput screening (experi- 
mental'”~*" or virtual***’). Although the efficacies 
of many of these inhibitors were supported by assay 
results, the modes of action are unknown for most 
of them. Also, the lack of structural information for 
inhibitors binding to MP" impedes the structure- 
based optimization of these inhibitors. To our 
knowledge, the crystal structure of MP"? bound by 
chloromethyl ketone (CMK) has so far been the 
only published structure for an inhibitor-bound 
MP". 

Aza-peptide epoxides (APEs) were synthesized 
as a new Class of inhibitors apparently specific for 
clan CD cysteine peptidases** (based on the 


weighting options of LEONORA,” yee the following 
parameters: Kinact=35(417) X10~ =18(+9) uM 
and K,,=96(+31) uM. reerea ae details 
are provided in Materials and Methods. 
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classification by Barrett and Rawlings™), including 
the legumains~° and the caspases.*” Each APE has 
an aza-peptide component, with an epoxide moiety 
attached to the carbonyl group of the P1 residue. 
The side-chain of the P1 residue predominantly 
determines the target-peptidase specificity of an 
APE. The substituent on the epoxide C2 atom also 
allows some tuning of both the inhibitory activity 
and specificity of APE towards a particular target 
peptidase. The aza-peptide component resembles a 
peptide, except that the C* atom of the P1 residue in 
the former is replaced by a nitrogen atom to form an 
aza-amino acid residue. This introduces trigonal 
planar geometry to the a-atom of the P1 residue and 
reduces the electrophilicity of the carbonyl C atom 
of the P1 residue, thereby making the carbonyl 
eroup of the P1 residue resistant to nucleophilic 
attack.*° It has been proposed that APEs inhibit 
their target peptidases irreversibly by a mechanism 
in which the catalytic Cys S” atom nucleophilically 
attacks one of the two epoxide carbon atoms (C2 or 
C3) of APE (Figure 1(a)).***°*” This results in the 
opening of the conformationally strained epoxide 
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ring, and the formation of a covalent bond between 
the Cys SY atom and the attacked APE atom. 


Results 


Inhibition of SARS-CoV MP"° by an APE 


We hypothesized that an APE possessing an aza- 
glutamine (AGIn) as the P1 residue to mimic the S1 
specificity of SARS-CoV MP" for Gln (Figure 1(a)) 
would irreversibly inhibit the peptidase. Accord- 
ingly, we synthesized Cbz-Leu-Phe-AGIln-(S,S)EP- 
COOEt and Cbz-Leu-Phe-AGIn-(R,R)EP-COOEt. 
Micromolar quantities of the S,S diastereomer 
strongly inhibited the cleavage of a peptidic 
substrate, manifesting itself as a pronounced 
slowing of the reaction velocity as the reaction 
progressed (Figure 1(b)). Under these conditions, 
the R,R diastereomer did not detectably inhibit 
MP". Analysis of the rates of inactivation (j,) at 
different concentrations of substrate and inhibitor 


Table 1. Parameters and statistics derived from X-ray diffraction data processing and structure refinement 


Unbound MP*° 

Data processing 
Wavelength (A). 1.116 
Resolution limit* (A) 50.00—2.08 (2.15—2.08) 
Space group C2 
Unit-cell constants 

a (A) 107.76 

b (A) 82.14 

c (A) 53.13 

a (°) 90.00 

B (°) 104.87 

y (°) 90.00 
Mosaicity (°) 1.00 


Total number of reflections 
Number of unique reflections 


50,960 (3872) 
26,196 (2482) 


Redundancy 1.9 (1.6) 
Completeness (%) 95.8 (91.0) 
Reym (%) 3.7 (22.5) 
I/o(J) 14.9 (3.4) 
Structure refinement. 

Resolution range (A) 40.00-2.08 
Rywork (%) 19.2 
Reree (%) 24.6 


MP"°: APE Mea): APE 
1.116 1.116 
50.00-1.88 (1.95-1.88) 50.00-2.30 (2.38-2.30) 
C2 P2122, 
106.71 65.35 
83.68 67.48 
52.87 167.32 
90.00 90.00 
105.66 90.00 
90.00 90.00 
0.61 0.23 


74,092 (6410) 
35,576 (3301) 


134,161 (10,755) 
33,305 (3147) 


2.0 (1.7) 4.0 (3.4) 
96.3 (90.0) 99.2 (95.5) 
3.3 (31.5) 5.2 (51.8) 
21.5 (2.1) 24.8 (2.0) 
40.00-1.88 40.00-2.30 

19.8 17.5 
24.2 24.9 


Number of non-hydrogen atoms per asymmetric unit (average B factor, A) 


Protein 2371 (56.03) 
APE Not applicable 
Solvent 150.5 (66.13) 


rms deviation from ideal geometry 


Bond lengths (A) 0.029 

Bond angles (°) 2.801 
Ramachandran plot 

Favored (%) 87.2 

Allowed (%) 10.6 

Generously allowed (%) 1.1 

Disallowed (%) 1 


* Numbers in the parentheses refer to the highest resolution bins. 


2371 (42.72) 
46 (45.28) 
264.5 (55.29) 


4690 (46.61) 
92 (65.29) 
279.5 (56.28) 


0.022 0.021 
2.342 1.944 
91:5 88.5 
6.8 29 
1.5 0.8 
0.4 0.8 


Rgym = Ene Zine, — Cn) \/ 2a Zne,i, Where Inia; and (Ipu1) are the ith observed intensity and average intensity of the reflection hkl, 


respectively. 


* Rwork = VIIFol —|F.II/Z1F,|, where |F,| and |F.| are the observed and calculated structure factor amplitudes of a particular reflection, 


respectively, and the summation is over 95% of the reflections in the specified resolution range. The remaining 5% of the reflections were 
randomly selected before the structure refinement and not included in the structure refinement. Rge~. was calculated over these 
reflections using the same equation as for Rain 
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Figure 2. SARS-CoV MP" in a C2 unit cell. The two 
protomers making up the dimeric peptidase are from two 
different asymmetric units (green and magenta) related 
by the crystallographic 2-fold symmetry along the b-axis 
(shown as arrows). 


indicated that APE inhibited MP" with a kinact/ Ki= 
1900(+400) M~'s—¢. 


Structure determination 


We have determined the crystal structures of 


SARS-CoV MP"’ in three forms: wild-type peptidase 
in the absence and presence of APE in the space 
group C2, and a variant of the peptidase with an Ala 
added to the N terminus of the wild-type sequence, 
Mi, bound by APE in the space group P22;21. 
The parameters and statistics derived from X-ray 
diffraction data processing and structure refine- 
ment are summarized in Table 1. For unbound MP" 
and the MP"°:APE complex in the space group C2, 
each asymmetric unit has only one protomer of the 
dimer. The two protomers of each dimer are related 
by the crystallographic 2-fold symmetry (Figure 2). 
All residues of the protomer (residues 1 to 306) were 
identified in the electron density maps. In the 
Ramachandran plot for the structure of unbound 
MP", Asp33, Ala46 and Glu47 are in the generously 
allowed regions, whereas Asn84, Tyr154 and [le286 
are in the disallowed regions. The Asp33 O* atom 
forms a hydrogen bond with the Tyr101 phenolic 
OH group (2.9 A). The Asn84 N°’ atom forms a 
hydrogen bond with the Glu178 carbonyl O atom 
(3.2 A), and possibly there are hydrogen bonds and 
van der Waals forces between the side-chains of 
Asn84 and Lys180 as well. Hydrophobic inter- 
actions occur between the side-chains of Thr285 
and Ile286 of opposite protomers at the dimer 
interface. The electron densities for the side-chains 
of Ala46, Glu47 and Tyr154 are not well defined. 
Similarly, in the Ramachandran plot for the 
structure of the MP":APE complex, Asp33, Asn84, 
Tyr154 and Asn277 are in the generously allowed 


Table 2. Root-mean-square differences (rmsd; in A) for superimpositions of the structures of SARS-CoV MP"® determined 


in different conditions 


Unbound MP"® (C2, pH 6.5) 
Unbound MP"® (C2, pH 6.5) 


MP'°: APE (C2, pH 6.5) 


Unbound MP" (P2,, pH 6.0) AP: 0.63 


Unbound MP" (P2,, pH 7.6) AP: 0.92 


Unbound MP"® (P2,, pH 8.0) AP: 0.92 


MP'°:CMK (P21, pH 6.0) AP: 0.58 


MP'°: APE (C2, pH 6.5) 


Mi aca: APE (P2;2)21, pH 6.5) 


PP: 0.24 PA: 0.35 
DD: 0.37 PB: 0.35 
DD: 0.65 

PA: 0.39 

— PB: 0.38 
DD: 0.74 

AP: 0.75 AA: 0.62 
BP: 0.37 BB: 0.39 
DD: 0.48 AB: 0.62 
BA: 0.34 

DD: 0.79 

AP: 0.99 AA: 0.76 
BP: 1.02 BB: 0.89 
DD: 1.04 AB: 0.82 
BA: 0.90 

DD: 0.83 

AP: 0.98 AA: 0.75 
BP: 1.01 BB: 0.85 
DD: 1.06 AB: 0.78 
BA: 0.84 

DD: 0.85 

AP: 0.68 AA: 0.58 
BP: 0.37 BB: 0.41 
DD: 0.51 AB: 0.57 
BA: 0.40 

DD: 0.70 


Superimpositions of structures were done for pairs of protomers (PA, AP, PB, BP, AA, BB and AB) and pairs of dimers (DD). The first 
letter corresponds to the part from the structure listed at the left of the Table superimposed and the second letter corresponds to the part 
from the structure listed along the top of the Table superimposed. P, protomer (applicable to C2 structures only); A, protomer A; B, 


protomer B; D, dimer. 
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region, whereas [le286 is in the disallowed region. 
The electron densities for Oe and Asn277 are 
not well defined. For the are igs complex in 
the space group P2)2,2), there is a dimer in each 
asymmetric unit. Only the residues 1A to 304A and 
residues 1B to 300B were identified in the electron 
density maps. In the Ramachandran plot for this 
structure, Asp33 and Asn84 of both protomers are 
in the generously allowed regions, whereas Tyr154 
and Ile286 of both protomers are in the disallowed 
regions. The electron density for Tyr154 is not well 
defined. Superimposition (see Materials and 
Methods) of protomers A and B yielded a root- 
mean-square difference (rmsd) of 0.29 A. Positional 
differences (up to 4.18 A) occur mainly among the N 
and C-terminal residues as well as those poorly 
defined residues on the flexible loops. The three 
structures are in close agreement (Table 2). With 
regard to the protomer orientation and protein fold, 
these structures, in general, are identical with the 
crystal structures of MP" previously determined in 
the space group P2,.'° 


Binding of APE to SARS-CoV MP"° 


Crystals of SARS-CoV MP and Mi',_,. were 
soaked in the solutions of Cbz-Leu-Phe-AGIn-(S,S) 
EP-COOEt, Cbz-Leu-Phe-AGIn-(R,R)EP-COOEt, 
and a racemic mixture of the S,S and R,R 


diastereomers (trans). Only the S,S diastereomer 
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showed up in the electron density maps. APE binds 
in the substrate-binding regions of MP*® (Figure 3(a) 
and (b)). As visualized in all three structures, the 
residues forming the substrate-binding regions of 
both protomers of the peptidase are in the 
catalytically competent conformation, similar to 
their counterparts in the structures of the main 
peptidases from other coronaviruses,'”*’ and to 
those in protomer A of the P2; structures of SARS- 
CoV MP?°.!S In the structure of unbound M?"’, the 
catalytic dyad has a distance of 3.7 A between 
the His41 N* atom and the Cys145 SY atom, and the 
Cys145 S* atom is coplanar with the atoms of 
the His41 imidazole ring. Superimposition of the 
structures of unbound MP" and the MP"°:APE 
complex shows that the binding of APE does not 
cause any major changes in the structure of the 
peptidase (Table 2). The Cys145 C’-C° bond under- 
goes a 95° rotation (x1: from —64.0° to +30.7°) 
accompanying the formation of a covalent bond 
with a distance of 2.01 A between the Cys145 S” 
atom of the peptidase and the epoxide C3 atom of 
APE (Figure 4(a) and (b)). The length of a C-S 
single bond is normally about 1.8 A. However, 
with an estimated overall coordinate error 
(based on maximum likelihood) of 0.12 A for 
the structure of the MP":APE complex, the dif- 
ference between the refined and _ expected 
distances (0.2 A) is not considered significant. 
This new covalent bond makes a torsion angle, 


Figure 3. Binding of APE (orange) in the substrate-binding regions of SARS-CoV MP". (a) Stereo view of the 
outstanding density in the F,—F, map for the structures of the MP"®:APE complex and protomer B of the Mia: APE 


pro 


complex. (b) The corresponding stereo view for protomer A of the M}),_,): APE complex. 
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Figure 5. Newman projections. (a) The equivalent to 6 of P1-AGIn, C(P2-Phe)-N(P1-AGIn)-N*(P1-AGIn)-C(P1-AGIn), 
(b) the equivalent to } of P1-AGIn, N(P1-AGIn)-N”*(P1-AGIn)-C(P1-AGIn)-C3(epoxide), (c) the torsion angle O=C(P1- 
AGIn)-C3(epoxide)-S*(Cys145), and (d) the torsion angle C(P1-AGIn)-C3(epoxide)-C2(epoxide)-C(ethyl ester carbonyl) 
of APE in the structure of the SARS-CoV MP"*: APE complex. (e) The equivalent to o of P1-AGIn, C(P2-Phe)-N(P1-AGIn)- 
N*(P1-AGIn)-C(P1-AGIn), (f) the equivalent to y of P1-AGlIn, N(P1-AGIn)-N’(P1-AGIn)-C(P1-AGIn)-C3(epoxide), 
(g) the torsion angles O=C(P1-AGIn)-C3(epoxide)-S*(Cys145) and (h) the torsion angle C(P1-AGIn)-C3(epoxide)- 
C2(epoxide)-C(ethyl ester carbonyl) of APE in protomers A and B of the Mera APE complex. 


O=C(P1-AGIn)-C3(epoxide)-S"(Cys145), of 65.4° 
(Figure 5(c)). The epoxide ring of APE opens, 
leaving a hydroxyl group on the C2 atom to form 
hydrogen bonds with the Asn142 O°! atom of the 


peptidase (2.9 A) and the P2-Phe carbonyl O atom 
of APE (3.2 A) (Figure 4(a) and (b)). The configu- 
rations of the C2 and C3 atoms are inverted from S,S 
to R,R. The conversion at the C3 atom is due to the 


Figure 4. (a) View from the floor of the S1 specificity pocket in all structures. (b) Schematic diagram for the interactions 
of APE with SARS-CoV MP" in the structures of the MP"°: APE complex and protomer B of the Miia(ay: APE complex. 
(c) The corresponding schematic diagram for protomer A of the M? Act): APE complex. Hydrogen bonds are shown as 
broken lines, with their distance(s) (in A) given in (b) and (c). In (b), the first value corresponds to the MP"’: APE complex 
and the second value corresponds to protomer B of the ME ay APE complex. The residues of the peptidase in contact 
with APE are shown as arcs. 
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formation of this new covalent bond, whereas that 
at the C2 atom is a result of the change in the 
ranking order of substituents (Figure 5(d)). Asn142 
and Cys145 of the peptidase flank a short loop 
constituting the oxyanion hole, with the amide 
hydrogen atoms of Gly143 and Cys145 as the 
hydrogen-bond donors to stabilize the carbonyl O 
atom in the scissile peptide group of the substrate. 
The conformation of this loop remains essentially 
unchanged upon the binding of APE. The P1-AGIn 
carbonyl O atom of APE is accommodated in the 
oxyanion hole of the peptidase, accepting hydrogen 
bonds from the amide hydrogen atoms of Gly143 
and Cys145 at distances of 2.7A and 3.44, 
respectively (Figure 4(a) and (b)). The P1-AGIn N® 
atom of APE remains sp*-hybridized and adopts 
trigonal planar geometry. For P1-AGIn, the equiva- 
lent to , C(P2-Phe)-N(P1-AGIn)-N%(P1-AGIn)- 
C(P1-AGIn), is —84.8° (Figure 5(a)), and the 
equivalent to y, N(P1-AGIn)-N*(P1-AGIn)-C(P1- 
AGIln)-C3(epoxide), is —7.3° (Figure 5(b)). With 
the error in atomic coordinates considered, the 
trigonal plane centered at the P1-AGIn N® atom is 
coplanar with that centered at the P1-AGIn carbonyl 
C atom, allowing the N® atom to reduce the 
electrophilicity of the carbonyl C atom by 
m-electron delocalization. The equivalent to y, of 
P1-AGIn, N-N*-C?-C”, is —77.4°. 

Compared to the structure of the MP*°:APE 
complex, the structure of the Mi): APE complex 
shows some differences in the geometry of binding. 
In the latter, the rotation of the Cys145 C*-C? bond 
of the peptidase reaches a more positive value of 71 
(protomer A, 47.5°; protomer B, 46.5°). The length 
of the covalent bond between the Cys145 S* atom of 
the peptidase and the epoxide C3 atom of APE is 
2.09 A in protomer A and 2.05 A in protomer B 
(Figure 4(a) to (c)). Note that the estimated overall 
coordinate error (based on maximum likelihood) 
for this structure is 0.17 A. O=C(P1-AGIln)- 
C3(epoxide)-S’(Cys145) makes a torsion angle of 
72.0° in protomer A and 85.5° in protomer B 
(Figure 5(g)). In both protomers, the configurations 
of the C2 and C3 atoms of APE are inverted from S,S 
to R,R (Figure 5(h)), and the P1-AGIn N” atom of 
APE remains sp*-hybridized and adopts trigonal 
planar geometry. For the P1-AGlIn of APE, the 
equivalent to , C(P2-Phe)-N(P1-AGIn)-N’(P1- 
AGIn)-C(P1-AGIn), is —83.7° in protomer A and 
—105.7° in protomer B (Figure 5(e)), and the 
equivalent to y, N(P1-AGIn)-N*(P1-AGIn)-C(P1- 
AGln)-C3(epoxide), is —12.3° in protomer A and 
— 8.8° in protomer B (Figure 5(f)). With the errors in 
atomic coordinates considered, the trigonal plane 
centered at the P1-AGlIn N® atom is roughly 
coplanar with that centered at the P1-AGIn carbonyl 
C atom in both protomers. The equivalent to x of 
P1-AGIn, N-N%-C?-C’, is —82.7° in protomer A and 
—119.5° in protomer B. 

The aza-peptide component of APE binds in the 
substrate-binding regions of both protomers of MP*® 
in a substrate-like manner. Consistent with the 
results of the secondary-structure studies for the 
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substrates of this peptidase,” the structures of the 
MP: APE and Mj',,:APE complexes show that 
the main chain of the aza-peptide component of 
APE interacts with the main chain of the residues 
164 to 166 of the peptidase through amide 
hydrogen-carbonyl oxygen hydrogen bonding in 
the manner of an anti-parallel B sheet in both 
protomers (Figure 4(a) to (c)). 


Specificity pockets 


The predominant S1 specificity of SARS-CoV 
MP"® for Gln is determined primarily by His163. In 
the structure of unbound MP'°, the His163 N® atom 
interacts with a chloride ion at a distance of 3.3 A 
and in the plane of the His163 imidazole ring, 
whereas in the structures of the MP*:APE and 
M\4y:4PE complexes, the chloride ion is dis- 
placed by the P1-AGIn side-chain amide group of 
APE with its O*' atom accepting a hydrogen bond 
from the His163 N* atom of the peptidase (2.6 A to 
2.8 A). Additional hydrogen bonds may be 
donated, though not of ideal geometry, by the 
P1-AGIn N* atom of APE to the Phe140 carbonyl O 
atom and the Glul66 O* atom of the peptidase 
(Figure 4(b) and (c)). In both protomers of all three 
structures, the Phe140 pheny] ring interacts with the 
His163 imidazole ring through 7-stacking (distance 
between the geometric centers of the rings: 3.7 to 
3.8 A), and the latter is properly oriented for its N°’ 
atom to accept a hydrogen bond from the Tyr161 
phenolic OH group (2.9 to 3.1 A). This hydrogen 
bond maintains the neutral tautomeric state of the 
His163 imidazole ring with its N°? atom protonated 
over a broad range of pH. This is crucial for the 
interaction of the His163 of the peptidase with the 
P1-Gln of the substrate.'°1°*” 

In each peptidase dimer, the integrity of the S1 
specificity pocket in one protomer requires the 
protonated amino group of Serl from the other 
protomer. This residue is at the tip of the N-finger 
(N-terminal residues 1 to 7) propagating between 
the domain II of its parent protomer and the 
domain II of the partner protomer. In the structures 
of unbound MP" and the MP":APE complex, the 
Serl amino group of each protomer forms an ionic 
interaction with the Glul66 side-chain carboxyl 
eroup of the other protomer (2.7 A). The Serl of 
each protomer also interacts with the Phe140 of the 
other protomer through amide hydrogen-—carbonyl 
oxygen hydrogen bonding (3.2 A). These three 
residues are thus held together to form the “floor” 
of the S1 specificity pocket (Figure 6(a)). In both 
protomers of the M!,_,:APE complex, however, 
Serl is N-terminally blocked by the additional Ala, 
and the ionic interaction between the Serl amino 
eroup and Glu166 side-chain carboxyl group is lost. 
This is also followed by the loss of the amide 
hydrogen-carbonyl oxygen hydrogen bond 
between Serl and Phe140 because Ser1 is no longer 
properly oriented. The additional Ala is disordered, 
leaving Serl unanchored and the floor of the S1 
specificity pocket partly disrupted (Figure 6(b)). 
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Figure 6. Floor of the S1 specificity pocket of SARS-CoV MP"’. Phe140, Glu166 and Glu290 of the first protomer (green) 
and the N-finger (N-terminal residues 1 to 7) of the second protomer (magenta) are shown. (a) In the structures of 
unbound MP?" and the MP"’:APE complex, Phe140 and Glu166 of the first protomer interact with Serl of the second 
protomer to form the floor of the S1 specificity pocket. (b) In the structure of the M®" A(t): APE complex, the disordered N- 
terminal Ala of the second protomer (not shown) orients Serl poorly for any interaction with Phe140 and Glu166 of the 
first protomer. The interaction between Glu290 of the first protomer and Arg4 of the second protomer, however, remains 
in both (a) and (b). Hydrogen bonds are shown as broken lines. 


Despite the disruption of the 51 pocket by the 
presence of a single additional N-terminal residue, 
the presence of a ten-residue affinity tag at the N 
terminus of MP™ reduced the specific activity of the 
peptidase by less than an order of magnitude 
(results not shown). 

Mr has greatest preference for Leu and lle at P2, 
followed by Phe, Val and Met in that order.” In the 
structures of the M?™:APE and Mia = APE 

complexes, the P2-Phe side-chain of APE fits snugly 
in the 52 specificity pocket of the peptidase, where 
the interactions are mainly hydrophobic. The 
P2-Phe phenyl ring of APE interacts with the 
His41 imidazole ring of the peptidase through 
m™-stacking (distance between the geometric centers 
of the two rings: 4.3 to 4.6 A). Superimposition of 
the structures of unbound MP" and the MP**:APE 
complex shows that, upon the binding of APE, the 
side-chain of Met49 undergoes a large confor- 
mational change, thereby opening the S2 specificity 
pocket for the P2-Phe side-chain of APE. Also, the 
side-chain of GIn189 is no longer disordered after its 
reorientation and formation of a hydrogen bond 
through its O*' atom with the P2-Phe amide 
hydrogen atom of APE (2.8 A to 3.4A). This 
appears to secure the P2-Phe of APE in the 82 
specificity pocket of the peptidase (Figures 3(a) and 
(b) and 4(a)). 

In the structures of the MP": APE and Mi",,_,): APE 
complexes, the P3-Leu side-chain of APE extends 
into the solvent and does not have any significant 
interactions with the peptidase (Figures 3(a) and (b) 
and 4(a)). This lack of interactions is consistent with 


the fact that no S3 specificity of MP" could be 
established.” 

The aza-peptide component of APE consists of 
only three residues, so its benzyloxycarbonyl (Cbz) 
eroup partly takes up the space for the P4 residue of 
a substrate. MP" has a shallow S4 specificity pocket 
that accommodates small side-chains (Ser, Thr, Val, 
Pro and Ala).’ In the M?"°:APE complex and 
protomer B of the Mj,,):APE complex, the 
binding of the Cbz group of APE does not use the 
54 specificity pocket of the peptidase. The position 
for the C* atom of a P4 residue of a substrate is 
occupied by the O2 atom in the Cbz group of APE. 
In this conformation, the benzyl group of APE 
makes contacts with Prol68 and with residues 190 
to 192 of the peptidase; the Cbz group is exposed to 
the solvent (Figures 3(a) and 4(b)). In contrast, in 
protomer A of the M{’,_,): APE complex, the benzyl 
eroup of APE squeezes into and thereby widens the 
54 specificity pocket of the peptidase, so that it is 
snugly accommodated in this enlarged pocket now 
formed by the residues 165 to 168, Phe185, Gln192 
and the main-chain atoms of Val186 (Figures 3(b) 
and 4(c)). 

The S1' specificity pocket of MP" is also shallow, 
accommodating only small side-chains (Ser, Ala, 
Gly, Asn and Cys).”” In the design of APEs, several 
different epoxide derivatives were attached to the 
aza-peptide component to modulate the inter- 
actions of APE with the S1’ specificity pocket of a 
clan CD cysteine peptidase.***°”” In the structures 
of the MP'°:APE and Mj,,:APE complexes, 
however, the S1’ specificity pocket of the peptidase 
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is apparently not used in the binding of APE. The 
epoxide C2 atom of APE sits close to the position for 
the C* atom of the P1’ residue of a substrate and the 
hydroxyl group on the C2 atom is exposed to the 
solvent. The ethyl ester group of APE lies against 
the “ceiling” of the active site lined by Leu27, Pro39, 
His41, and the peptide group between His41 and 
Val42. 


Dimer interface 


The structures reported here show similar 
features at their dimer interfaces. The solvent- 
accessible surface area (per protomer) buried 
upon dimerization of SARS-CoV MP" is from 
1250 A* to 1260 A*. In the structures of unbound 
MP’ and the MP'°:APE complex, the crystallo- 
eraphic 2-fold axis passes through the dimer 
interface and brings the opposite interacting 
residue-pairs into exact 2-fold symmetry (Figure 2); 
whereas in the structure of the M\,,:APE 
complex, even in the absence of a crystallographic 
2-fold axis, the dimer interface still exhibits 
approximate 2-fold symmetry. In all three struc- 
tures, the dimer interface concentrates on one face 
of each protomer: that containing the residues on 
the N-finger, and domains II and HI. The majority of 
the interactions occur between the residues on the 
N-finger and domain II. Two ionic interactions, one 
between the Serl amino group and the Glu166 side- 
chain carboxyl group of opposite protomers, and 
one between the side-chains of the Arg4 and Glu290 
of opposite protomers, are observed in the struc- 
tures of unbound MP" and the Mi,_,):APE 
complex (Figure 6(a)). However, in the structure of 
the Mi",_,): APE complex, only the latter ion pair is 
observe (Figure 6(b)). In contrast to the N-finger 
immobilized at the dimer interface, the C-terminal 
loop (residues 301 to 306) is highly mobile as it is 
exposed to the solvent and anchored by only three 
or four solvent-exposed hydrogen bonds to some 
residues along the rim of the dimer interface. 
Interestingly, in the structures of unbound MP" 
and the MP"®:APE complex, the C-terminal loop of 
one protomer extends towards the S1 specificity 
pocket of the other protomer, whereas in the 
structure of the M)’,_,:APE complex, the 
C-terminal loop of protomer A (the only one 
visualized in this structure) propagates away from 
the S1 specificity pocket of protomer B. 


Discussion 


The kinetic data and crystal structures reported 
herein indicate that APEs have excellent potential 
as inhibitors of SARS-CoV MP" and are worthy of 
further evaluation in the development of lead 
compounds for anti-SARS agents. Thus, the Kinact/ 
K; of MP*° for Cbz-Leu-Phe-AGIn-(S,S)EP-COOEt is 
similar in magnitude to that of the first generation 
APEs produced to inhibit other cysteine pepti- 
dases.***°*” Optimization of the latter has yielded 
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inhibitors of caspases with kijat/K, well over 
10°M*s 1.” 

The excellent specificities of APEs for clan CD 
cysteine peptidases*”*°~” suggest that APEs have 
better potential as inhibitors of MP" than do 
chloromethyl ketones (CMKs), the first class of 
potential inhibitors proposed on a structural basis. 
The structure of the MP"°’:CMK complex previously 
determined in the space group P2, shows different 
and unexpected modes of binding for CMK to the 
two protomers of the peptidase.'? CMKs are highly 
active alkylating agents and react with good 
nucleophiles, such as hydroxyl and thiol groups. 
They therefore inhibit serine peptidases as well as 
cysteine peptidases.~* A recent study showed that 
CMKs efficiently inhibit some clan CA cysteine 
peptidases, such as papain and the cathepsins.” 
This casts doubt on the utility of CMKs as specific 
inhibitors for MP*°. In contrast, the structures of 
unbound MP" and the MP**:APE complex show 
that the aza-peptide component of APE binds to the 
peptidase in a substrate-like manner. The main- 
chain of the aza-peptide component of APE forms 
amide hydrogen-carbonyl oxygen hydrogen bonds 
with the main chain of the residues 164 to 166 of the 
peptidase in the manner of an anti-parallel B-sheet. 
The P1 and P2 side-chains of APE occupy the 51 and 
52 specificity pockets of the peptidase, respectively. 
Based on the definitions for the binding of 
epoxysuccinyl peptides to clan CA cysteine pepti- 
dases,*° this corresponds to the S and S’ binding 
mode, with inclination to the S binding mode 
because the pre-cleavage portion of the substrate- 
binding region of the peptidase makes the major 
contribution to MP"°:APE interactions. 

The structures of the MP"°:APE and Mi',_,): APE 
complexes substantiate the mechanism by which 
APEs have been proposed to irreversibly inhibit 
their target peptidases (Figure 1(a)). Nucleophilic 
attack at the epoxide C3, rather than the C2, atom of 
APE, by the Cys145 S” atom is consistent with the 
expected transition-state geometry for proteolysis 
catalyzed by MP". In caspase-3, the epoxide C3 
atom is attacked (M. Grutter, unpublished results), 
whereas in caspase-1, the C2 atom is attacked 
(R. Rubin, unpublished results). In the case of 
epoxysuccinyl peptides binding to clan CA cysteine 
peptidases, the position of attack depends on the 
orientation of the epoxysuccinyl peptide in the 
substrate-binding region.** E-64 binds to papain in 
the S binding mode, similar to the mode of APE 
binding to MP**. In the papain:E-64 complex, 
however, the epoxide C2 atom is the one attacked.”’ 
Nucleophilic attack at the epoxide C3 atom is 
observed in the S’ binding mode, as exemplified 
by CA-074 binding to cathepsin B.** 

The kinetic data and crystal structures indicate 
that MP" reacts only with the S,S diastereomer of 
the APE and not its R,R diastereomer. Interestingly, 
the order of inhibitory activities for APEs towards 
most clan CD cysteine peptidases is $5,5>R,R> 
trans>cis (the racemic mixture of the S,R and R,S 
diastereomers).***°*’ Based on the MP*°:APE 
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Figure 7. Models for each of the four diastereomers of Cbz-Leu-Phe-AGIn-EP-COOEt binding to SARS-CoV MP*® 
before the nucleophilic attack by the Cys145 S* atom of the peptidase. (a) 2S, 3S, (b) 2R, 3R, (c) 2S, 3R, (d) 2R, 3S. 


structures, we built models of each of the four 
diastereomers of APE at the active site of MP*° to 
explain the results of our trials (Figure 7(a) to (d)). 
These models show that, for APE to be accommo- 
dated in the substrate-binding regions of the 


peptidase, the epoxide C3 atom of APE must be in 
the S configuration, otherwise the epoxide moiety 
sterically clashes with the “back-wall” of the active 
site of the peptidase and with the aza-peptide 
component of APE itself (Figure 7(b) and (c)). In the 
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S configuration, the epoxide C3 atom of APE is also 
in better geometry with respect to the Cys145 S” 
atom for the nucleophilic attack. The epoxide C2 
atom of APE must be in the S configuration as well 
to allow the interactions between the epoxide 
moiety and the active site of the peptidase, 
otherwise the epoxide moiety sterically clashes 
with the loop constituting the oxyanion hole of the 
peptidase (Figure 7(d)). The model with the S,S 
diastereomer of APE also indicates that the distance 
between the His41 N* atom of the peptidase and 
the epoxide O atom is 4 A to 5 A, and that these two 
atoms are not well aligned for proton transfer 
(Figure 7(a)). This suggests that the opening of the 
epoxide ring likely involves two steps separated by 
a conformational rearrangement of APE: (1) the 
protonation of the epoxide O atom and (2) the 
nucleophilic attack at the epoxide C3 atom. From 
the model, it is not possible to determine the order 
in which these steps occur. However, it would be 
energetically more favorable for protonation to 
be the first step. On the other hand, MP"° may be 
sufficiently flexible in solution to allow the align- 
ment of the His41 N** atom and the epoxide O atom 
for proton transfer. In such a scenario, protonation 
and nucleophilic attack could occur in a concerted 
manner, enabling the epoxide ring to open in a 
single step. Much of the mechanism for the 
inhibition of SARS-CoV MP"® by APE remains to 
be elucidated. Rigorous treatment of this issue 
using methodologies in organic chemistry will be 
required. 

All three structures successfully visualize the 
N-fingers of both protomers of MP"’. Despite the 
appreciable participation by the N-fingers in dimer 
interactions, it was shown experimentally that the 
deletion of the N-fingers inactivates the peptidase 
but has little effect on its dimerization properties. 
Molecular dynamic simulations suggested that the 
main role of the N-finger is one directing the 
peptidase to dimerize at an orientation facilitating 
the formation of substrate-binding regions in the 
catalytically competent conformation. This very 
likely relies on the two ionic interactions observed 
in the structures of unbound MP" and the 
MP":APE complex, one between the protonated 
amino group of Serl and the Glul66 side-chain 
carboxyl group of opposite protomers, and one 
between the side-chains of the Arg4 and Glu290 
of opposite protomers. In the structure of the 
Mi41:SPE complex, the former ion-pair does 
not exist but both substrate-binding regions are still 
capable of accommodating APE in a manner similar 
to that exhibited by the structure of the MP*’:APE 
complex. This suggests that the former ionic 
interaction, possibly because it is relatively acces- 
sible to the solvent and therefore weaker, is of less 
importance. 

The crystal structure of MP" previously deter- 
mined in the space group P2; at pH 6.0 showed the 
collapse of the active site and 51 specificity pocket 
of one of the protomers, whereas the P2, structures 
at pH 7.6 and 8.0 showed the recovery of the 
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collapsed parts. Based on this trend, a pH-triggered 
switch for the catalytic activity of the peptidase was 
proposed.'” In our study, all crystals were grown at 
pH 6.5 and the resulting structures show that the 
substrate-binding regions of both protomers are in 
the catalytically competent conformation. This 
suggests the possibility of an alternative or 
additional mechanism underlying the pH-depen- 
dence of the activity of the peptidase, especially at 
pH 6.5 or above. The change in protonation/ 
deprotonation state of the catalytic dyad with pH 
is one of the possible second mechanisms. Insights 
into this possibility could be provided by the direct 
determinations of the pK, values of the catalytic 
dyad by nuclear magnetic resonance (NMR) spec- 
troscopy, as exemplified by the similar studies of the 
catalytic triad of a-lytic peptidase.” 


Materials and Methods 


Preparation of the protein and APEs 


SARS-CoV ME eas was cloned, overexpressed and 
purified as described.*! A clone expressing MP*° was 
generated using oligonucleotide-directed evolution to 
delete the codon corresponding to the N-terminal Ala of 
Mi',4) Using this clone, MP*° was overexpressed and 
purified essentially as described for M'";_,.. Cbz-Leu- 
Phe-AGIn-(S,S)EP-COOEt and Cbz-Leu-Phe-AGIn- 
(R,R)JEP-COOEt were synthesized using the methods 
established to synthesize other APEs***°*” with minor 
modifications. 


Inhibition studies 


Enzymatic activity was measured by following the 
increase in fluorescence due to the cleavage of a 
fluorogenic peptide: Abz-Ser-Val-Thr-Leu-Gln-Ser-Gly- 
(NO2)Tyr-Arg, where Abz is aminobenzoate and 
(NO2)Tyr is nitrotyrosine. Fluorescence was measured 
using a Cary Eclipse Fluorescence spectrophotometer 
(Varian Canada, Mississauga, Ontario, Canada) equipped 
with a circulating water-bath. Experiments were per- 
formed using a 100 pl quartz cuvette. The standard assay 
contained 25nM MP", 20 mM Bis-Tris (pH 7.0), 2mM 
DTT, and was performed at 37.0(+0.1) °C. The reaction 
was monitored using an excitation wavelength of 320 nm 
(5 nm bandpass) and an emission wavelength of 420 nm 
(10 nm bandpass). Initial velocities were determined from 
a least-squares analysis of the linear portion of the 
progress curves (at least 1 min) using Excel 2003 
(Microsoft, Redmond, WA). All rates were corrected for 
the inner filter effect using an empirical correction.” 

In inhibition studies, the concentration of APE was 
varied from O0uM to 10 uM and the concentration of 
peptidic substrate was varied from 16 uM to 100 uM. 
These substrate and inhibitor concentrations were 
dictated by solubility limitations and the observed rates 
of inhibition. The rate of inactivation at each concen- 
tration of substrate and inhibitor, j,, was determined by 
fitting equation (1)*° to the corresponding progress curve 
using SCIENTIST version 2.01 (Micromath Scientific 
Software, Salt Lake City, UT). The parameters of 
inactivation, Kinace and K;, were evaluated by fitting 
equation (2) to the j, obtained at each concentration of 5 
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and I,°° using the least-squares and dynamic weighting 
options of LEONORA.” 


Crystallization, crystal soaking and cryo-protection 


Before crystallization, both SARS-CoV MP" and 
ME cs were dialyzed against 20mM NaCl, 20 mM 
Tris-HCl (pH 7.5), and concentrated to 10 mg/ml. All 
crystals were grown at ambient temperature by the 
hanging-drop, vapor-diffusion method. For the C2 
crystals, the reservoir solution contained 50mM 
ammonium acetate, 5% (w/v) polyethylene glycol (M, 
10,000), 3% ethylene glycol, 3% dimethyl sulfoxide, 1 mM 
dithiothreitol and 0.1mM Mes (pH 6.5). The drop 
contained equal amounts of the MP" solution and the 
reservoir solution. Block-shaped crystals grew in two to 
three days to a size of about 0.1 mm X0.1 mm X0.1 mm. 
For the P2,2,2, crystals, the reservoir solution had 
essentially the same composition as that for the C2 
crystals, except the replacement of 5% polyethylene 
glycol (M, 10,000) and 3% dimethyl sulfoxide by 6% 
polyethylene glycol (M, 8000). The drop contained equal 
amounts of the NV nas solution and the reservoir 
solution. Needle-shaped crystals grew in three to five 
days to a size of about 0.05mmX0.05 mm X0.5 mm. 
Crystals of good quality were selected and soaked 
overnight in drops with the same compositions as their 
reservoir solutions plus the APE chosen for this study at 
3 mM. Cryo-protectants had essentially the same compo- 
sitions as reservoir solutions, except for the inclusion of 
25% (v/v) glycerol in the case of the C2 crystals and the 
increase of ethylene glycol to 25% in the case of the 
P2,2,2, crystals. Crystals were briefly soaked and then 
immediately frozen in liquid nitrogen for storage and 
shipment to the synchrotron beamline. 


Data collection and processing, structure solution, 
refinement and analysis 


The X-ray diffraction data from all crystals were 
collected at the synchrotron Beamline 8.3.1 (equipped 
with an ADSC-Q210 CCD detector) at the Advanced 
Light Source in the Lawrence Berkeley National Labora- 
tory. All data sets were indexed, scaled and merged using 
DENZO and SCALEPACK.”® Structure solution and 
refinement were carried out in CCP4.°"*? All structures 
were solved by the molecular replacement method,”' 
using the structure of unbound SARS-CoV MP" at pH 8.0 
in the space group P2; (PDB accession code 1UK2)’” as 
the search model for the structure of unbound MP"® in the 
space group C2, and the structure of unbound MP" in 
the space group C2 as the search model for the structure 
of the MP":APE complex in the space group C2 and 
the structure of the Mica: APE complex in the space 
group P2;2,2;. The structures of unbound MP" and the 
MEA :APE complex were solved using AMoRe,* and 
the structure of the Mi ay: APE complex was solved 
using MOLREP.** APE was located as outstanding 
electron densities in the substrate-binding region of the 
peptidase in both the F, — F, (contoured at 30 and 40) and 
2F,—F, (contoured at 1o) maps for each structure. All 
structures were iteratively refined using REFMAC™ and 
adjusted using XtalView/Xfit.*° The stereochemical 
qualities of the final structures were assessed using 
PROCHECK.”*° Graphical representations of the struc- 
tures were prepared using PyMOL?7. Superimpositions of 
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structures were carried out using ALIGN,*”*® based on 
the main-chain atoms (amide N, C*, and carbonyl C and 
O). The surface areas of structures were calculated using 
NACCESS.£ APE-peptidase interactions and dimer 
interactions were analyzed using LIGPLOT and DIM- 
PLOT,” respectively. 


Protein Data Bank accession codes 


The atomic coordinates and structure factors of all 
structures have been deposited in the RCSB Protein Data 
Bank. The accession code is 2A5A for the structure of 
unbound SARS-CoV MP*®, 2A5I for the structure of the 
MP":APE complex and 2A5K for the structure of the 
ME AVE complex. 
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